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^ ! Abstract 

, A new general decomposition theory inspired from modular graph decomposition 



is presented. This helps unifying modular decomposition on different structures, 
including (but not restricted to) graphs. Moreover, even in the case of graphs, 
this new notion called homogeneous modules not only captures the classical graph 
modules but also allows to handle 2— connected components, star-cutsets, and other 
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^ ' vertex subsets. 

The main result is that most of the nice algorithmic tools developed for modular 
PsJ ■ decomposition of graphs still apply efficiently on our generalisation of modules. 

^ , Besides, when an essential axiom is satisfied, almost all the important properties can 

be retrieved. For this case, an algorithm given by Ehrenfeucht, Gabow, McConnell 
and Sullivan [15] is generalised and yields a very efficient solution to the associated 
decomposition problem. 
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1 Introduction 

■ 

•w^ - 

X 

■ Modular decomposition has arisen in different contexts as a very natural opera- 

tion on many discrete structures such as graphs, directed graphs, 2-structures, 
automata, boolean functions, hypergraphs, and matroids. In graph theory, 
modular decomposition plays a central role. Not only modular graph decom- 
position yields a framework for the computation of all transitive orientations 
of a given comparability graph [18,20,27], but it also highly relates to com- 
mon intervals of a set of permutations [2,6,30] and therefore has applications 
in bioinformatics. Besides, many graph classes such as cographs, P4-sparse or 
P4-tidy graphs are characterised by properties of their modules (see e.g. [4]). 
It is also worth noticing that well-known NP-hard problems such as colouring 
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Fig. 1. Illustration of classical modular graph decomposition, i. In this undirected 
graph, vertex c is a splitter of {a,b} (not linked the same way), whereas vertex b 
is not a splitter of {a, c}. Vertex set {d,e, f, g}, as well as any of its subsets, is a 
module of this undirected graph, ii. Modular decomposition tree of the graph. 

can be solved in polynomial, and often linear, time when the graph is "suffi- 
ciently" decomposable [29] using some application of the divide and conquer 
paradigm. Finally, the decomposition is useful for graph drawing [32] , compact 
encoding (e.g. with cographs [10] and P4— sparse graphs [25]), and precom- 
puting for graph problems including recognition, decision, and combinatorics 
optimisations (see [29] or [4] for a survey). A central point of this theory re- 
lies on the decomposition theorem which presents a tree, so-called modular 
decomposition tree, as compact encoding of the family of modules of a graph. 
Then, computing this tree efficiently given the graph has been an important 
challenge of the past three decades [6,7,8,9,11,13,15,21,23,27,29,30]. 

On the other hand, several combinatorial algorithms are based on partition 
refinement techniques [22,23,31]. Many graph algorithms make intensive use 
of vertex splitting, the action of splitting parts according to the neighbourhood 
of a vertex. For instance, all known linear-time modular decomposition algo- 
rithms on graphs use this technique [6,8,11,13,15,22,23,27]. In bioinformatics 
also, the distinction of a set by an element, so-called splitter, seems to play an 
important role, e.g. in the efficient computation of the set of common intervals 
of two permutations [6,33]. 

An abstract notion of splitter is studied here and a formalism based on the 
concept of homogeneity is proposed. The resulting structures will be referred 
to as homogeneous relations. Our aim is a better understanding of existing 
modular decomposition algorithms by characterising the algebraic properties 
on which they rely. As a natural consequence, the new formalism unifies modu- 
lar decomposition on graphs and on their common generalisations to directed 
graphs [26] and to 2— structures [16]. Of course, the theory still applies on 
structures beyond the previous ones. Moreover, even in the case of graphs, 
this new notion called homogeneous modules not only captures the classical 
graph modules but also allows to handle other vertex subsets, e.g. those similar 
to 2— connected components, or to star-cutsets. 
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Our main result is that most of the nice algorithmic tools developed to com- 
pute the modular decomposition tree of a graph still apply efficiently in the 
general theory. For graph modules, to design efficient algorithms there actu- 
ally are three main approaches, distinguishable by the use of properties of: 
the set of maximal modules excluding a vertex [15], a factoring permuta- 
tion [6,8,21,23], or the visit order of some pecuhar graph search such as the 
so-called LexBFS lexicographic breadth-first search [5,14,24]. Because of its 
specificity due to exotic graph searches, the use of the third approach in the 
new theory is forfeit. Still, we extend the two first approaches, and retrieve 
most of the common efficient computations. 

However, as a consequence of their broadness, no obvious decomposition theo- 
rem, to our knowledge, is available for arbitrary homogeneous relations, hence 
no homogeneous modular decomposition tree necessarily is guaranteed. In- 
deed, though the homogeneous modules inherit many interesting properties 
from graph modules, they do not necessarily satisfy the following essential 
one. One can shrink a whole graph module M into one single vertex m e M: 
if some vertex of M distinguishes two exterior vertices, then so does every 
vertex of M and so does m. Let us denote the property by the name of mod- 
ular quotient. It actually is the basis of many divide-and-conqucr paradigms 
derived from the modular graph decomposition framework, such as the com- 
putation of weighted maximal stable or clique set, and graph colouring [19,28]. 
This naturally motivates us to study homogeneous relations fulfilling the mod- 
ular quotient property, hereafter denoted by good homogeneous relations. As 
expected, almost all important properties of modular graph decomposition, in- 
cluding the decomposition theorem, still hold for the latter relations. Eventu- 
ally, we generalise an algorithm given by Ehrcnfeucht et al. [15] to an 
algorithm computing the decomposition tree of a given good homogeneous 
relation on X. 

The paper is structured as follows. First the new combinatorial decomposi- 
tion theory is detailed in Sections 2 and 3. Section 4 investigates the general 
algorithmic framework on arbitrary homogeneous relations. The subsequent 
Section 5 is devoted to good homogeneous relations. Finally, we close the paper 
with noteworthy outcomes. 



2 Homogeneity, an abstraction of Adjacency 

Throughout this section X is a finite set, and denotes the family of all 

subsets of X. A diverse triple is {x,y,z) C X^ with x ^ y and x ^ z. This 
will be denoted by {x\yz) instead of (x, y, z) since the first element plays a 
particular role. Let H he a, relation over the diverse triples of X. Given a; G X, 
we define as the binary relation on X \ {x} such that Hx{y, z) H[x\yz). 
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Fig. 2. The standard homogeneous relation H of this directed graph satisfies 
H{x\uiVi) for all i, and -^H{x\uiVj) for all i / j. 



Definition 1 (Homogeneous Relation) H is a homogeneous relation on 
X if, for all x & X , is an equivalence relation on X\ {x} (i.e. it fulfils the 
symmetry, reflexivity and transitivity properties). Equivalently, such a relation 
can he seen as a mapping from each x & X to a partition of X \ {x}, namely 
the equivalence classes of H^. 

Definition 2 (Homogeneous Module) Let H be a homogeneous relation 
on X. A subset M C X is a homogeneous module of H if 

Vm, m' e M, yxeX\M, H{x\mm'). 

Remark: From the definition it is obvious that, given a homogeneous module 
M, if -iH[x\mm') for some m,m' & M then x e M. 

If -iH{x\mm') we say that x distinguishes m from m', or x is a splitter of 
{m, m'}. A homogeneous module M is trivial if \M\ < 1 or M ~ X. The family 
of homogeneous modules of H is denoted by Mj/, and M when no confusion 
occurs. H is modular prime if M,h is reduced to the trivial homogeneous 
modules. For convenience, such a relation is also called prime when it clearly 
appears in the context that modules are involved. Homogeneity and distinction 
can be applied to graphs. Indeed, there is a natural homogeneous relation 
associated to graphs as follow. 

Definition 3 (Standeird Homogeneous Relation) The standard homoge- 
neous relation H{G) of a directed graph G — {X, A) is defined such that, for 
all x,u,v e X , H{G){x\uv) is true if and only if the two following conditions 
hold: 

1. either both u and v or none of them are in-neighbours of x, and 

2. either both u and v or none of them are out-neighbours of x. 
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Roughly, H{x\uv) tells if x "sees" u and v the same way. Of course the above 
definition also holds for undirected graphs, tournaments, oriented graphs, and 
can also be extended to 2— structures (which roughly are edge-coloured com- 
plete directed graphs G = {X,X^), see e.g. [16] for further information). It 
follows straight from definition that 

Proposition 1 Let G be a graph, resp. tournament, oriented graph, directed 

graph, 2— structure. Homogeneous modules of its standard homogeneous rela- 
tion H{G) are modules of G in the usual sense [16,18,29]. 

Standard homogeneous relations are closely related to the notion of adjacency 
in graph theory. Notice that there are other homogeneous relations bound to 
a graph or to a 2— structure (e.g. in Section 6). Let us now give some first 
structural properties of homogeneous relations. Given A <^ X one can define 
the induced relation H[A\ as H restricted to diverse triples of A^. li A is a, 
homogeneous module we have the following nice property: 

Proposition 2 (Restriction) Let H be a homogeneous relation, M a homo- 
geneous module of H, and N C M. Then, N e 'Mh[m] N & 'Mh- 

Proof: That a homogeneous module of i/ is a homogeneous module of H[M] is 
straight from definition. Conversely, if C M is not a homogeneous module of 
H, then there is a splitter s E X\N such that ^x,y G N, -^H^slxy). However, 
s cannot belong to X\M since this would imply s is a sphtter w.r.t. H of M. 
Therefore, s E M \ N, and is a sphtter w.r.t. H[M] of N. Hence, N is not a 
homogeneous module of H[M]. □ 

2.1 Lattice Structure 

Let H be an arbitrary homogeneous relation over a finite set X. Let M denote 
the family of its homogeneous modules. Two sets A and B overlap if An B, 
A\B and B\A all are non-empty. It is denoted by A^B. 

Proposition 3 VA, B e M, if A^B, then (A n 5) e M and {A\J B) E M. 

Proof: the fact that An B is a homogeneous module is obvious. We use the 
transitivity of for all x ^ ^ U S to prove (A U S) e M. □ 

Proposition 4 If M denotes the family of homogeneous modules of a homo- 
geneous relation, and M' = M U {0}, then (M', C) is a lattice. 

Proof: Since G M', and thanks to Proposition 3, the intersection of two 
members A and B belonging to M' belongs to M'. It is the infimum of A and 
B, since any member of M' that is a subset of both A and 5 is a subset of 
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An B. Let K be the family of all members of M' containing both A and B. 
It is non-empty for X is a member. Since M' is closed under intersection, 3\f 
admits a unique smallest member (w.r.t. inclusion), which is the intersection 
of all its members, and is the supremum of A and B. □ 

This lattice is a sublattice of the boolean lattice (hypercube) on X. Moreover, 
if we consider A eM such that \A\ > 1, and M(>1) = {M G Mh and M D A}, 
then (M(^), C) is a distributive lattice. 

2.2 Homogeneous Modules as Roots of a Submodular Function 

Submodular functions are combinatorial objects with powerful potential (see 
e.g. [17]). Theorem 1 below enables the application of this theory to homoge- 
neous relations: the homogenenous modules of any such relation coincide with 
the roots of a function which satisfies the submodular inequality on intersect- 
ing subsets. 

Definition 4 A set function ^ : J'(X) R submodular if, for all sets 
A,B e 9{X), ii{A) + ii{B) > ii{A UB)+ ii{A n B) (see e.g. [17]). 

Theorem 1 Let H he a homogeneous relation on X. Let s{A) he the function 

counting the number of splitters of a non-empty subset A C X. Then, s follows 
the submodular inequality on intersecting subsets: 

s{A) + s{B) > s(A UB) + s(A n B) for allAr\B^^. 

Proof: \i A <Z B or B <Z A, the inequality is trivial. If Ago 5 then A 7^ and 
B $. Let §A denote the set of splitters of A. If {Xi, . . . , Xk} is a partition of 
X, we note X = {Xi, Xk}. Obviously, S^ns = {§>AnB \ B, S^ns n B}. As 
Sa n A = 0, the partition §aub = {^aub \ ^a, ^aub H S^} can be reduced to 
§AuB = {^AuB \ §A, §A \ U B)}. Similarly, §b = {§b \ §AnB, §AnB \ B} . 
Finally, §a = {Sa \ B, {§a ^ B)\ Sahb, (§a H i?) fl Sadb} can be reduced to 
$A = {Sa \ {A U B), {§A n B) \ SahB: Bahb n B}. Hence, 

|§a| + |Sb| — |§aub| — |§AnB| = n S) \ SapibI + |Sb \ SapibI — |Saub \ 

To achieve proving the theorem, we prove that Saus\Sa Q ^bX^adb- Indeed, 
let s G Saub\Sa- Then, s ^ AUB and H{s\xy) for all x,y E A. Now, suppose 
that s ^ Sb- Since s does not belong to B, we deduce H{s\xy) for all x,y & B. 
Furthermore, as A and B overlap and thanks to the transitivity of H, we 
deduce H{s\xy) for all x,y G AU B and s ^ A U B, which is by definition 
s ^ §aub- Contradiction. Finally, supposing s G §AnB would imply s G S^- ^ 

In [33] a (restricted) version of this theorem is proved, and this submodularity 
property is used to propose a very nice algorithm which computes the set 
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Fig. 3. The atoms Ai, . . . ,Aq of the overlap class C = {Ci, C2, C3, C4}. 

of common intervals of a set of permutations. This approach was generalised 
for modules of standard homogeneous relations of undirected graphs in [6] . It 
would be interesting to consider this idea on arbitrary homogeneous relations. 

2.3 Strong Homogeneous Modules and Primality 

In an arbitrary family 3^ of subsets of X, a member A E 3^ is strong if it does 
not overlap any other member B eJ. Those which arc not strong arc weak. If 
they belong to the family, X and the singletons {x} {x G X) form the trivial 
strong members of 3^. Otherwise we extend 3^ with the trivial strong members. 

The set inclusion orders the strong members of 3^ into a tree, hereafter denoted 
by the generalised decomposition tree of 3^. This could be seen as a quick proof 
that, in 3^, there are at most 2\X\ — 1 strong members, and at most |X| — 2 
non-trivial ones, since the tree has \X\ leaves and no degree 2 internal nodes, 
except for possibly the root. When 5" is weakly partitive (see definition in 
Section 3), this tree plays an important role since it is an exact coding in 
space of the possibly 2'"^' members of the family. It is then called the 
decomposition tree of 5". 

The parent of a (possibly weak) member M G 3" is the smallest strong member 
Mp properly containing M, and M is said to be a child of Mp. For instance, 
if M is strong then Mp is its parent in the generalised decomposition tree. A 
strong member is prime if all its children are strong, and brittle otherwise. 

An overlap class of 5" is an equivalence class of the transitive closure of the 
overlap relation on 3". Such a class is trivial if it contains only one member 
A e 3^. Then A is by definition a strong member of 3^. The support of an 
overlap class C — {Ci, . . . , C^} is defined as S{C) = Ci U • • • U C^. An atom 
of the overlap class C is a maximal subset of S{C) that does not overlap any 
Ci {1 < i < k) (an illustration is given in Fig. 3). Notice that the atoms form 
a partition of S{C). Besides, an atom of an overlap class belongs to the class 
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if and only if this class is trivial. Furthermore, the support, resp. an atom, of 
an overlap class belongs to the family 5" if and only if it is a strong member 
of 5". Of course, a support, resp. an atom, does not necessarily belong to 9^. 
However, in a weakly partitive family (see Section 3), all atoms and supports 
of overlap classes will belong by definition of partitivity to (F, hence are strong 
members of the family. It is an elementary result of finite set theory that 

Proposition 5 The following holds for any family 3^ of subsets of a finite set 
X satisfying the closure under union of overlapping members. 

1. A C. X is a prime strong member of 3^ if and only if {A} is a trivial overlap 
class of 3^. 

2. A C X is a brittle strong member of 3^ if and only if it is the support of some 
non-trivial overlap class C a of 3^. In this case, weak children of A coincide with 
members of Ca- 

Of course we apply all these notions to the family of homogeneous modules 
of a homogeneous relation H. Let Z{x, y) be the largest homogeneous module 
of H containing x but not y. Z[x,y) is well defined since it is the union of 
all homogeneous modules containing x but not y, which is a homogeneous 
module thanks to Proposition 3. Moreover, Z{x, y) is not empty because {x} 
is a member. Let 2{H) be the family 

Z{H) = {Z{x,y) I x,yeX A x^y}. 

Notice that Z{H) is not necessarily closed under union of overlapping mem- 
bers. An example of such Z[H) is as follows. If X = (a, 6, c}, H{a\bc), H{b\ac), 
and H{c\ab), then {a, 6} e Z{H), {a, c} G Z{H), however X ^ Z{H). 

Theorem 2 All support and atoms of Z{H) that are homogeneous modules of 
H are strong homogeneous modules. A non-trivial strong homogenous module 
of H is either the support or an atom of some overlap class of Z{H). 

Proof: Let us prove the first claim of the theorem. 

(1) The support of an overlap class of Z{H) is a homogeneous module, since 
the family of homogenous modules is closed under the union of overlap- 
ping members (Proposition 3). If the support 5" of a given overlap class 
C is overlapped by another homogenous module, then it is overlapped by 
a homogenous module A ^ Z{H). Let x be an element of A \ 5* and y 
an element oi S \ A. Z{x,y) contains A but not {y} and thus overlaps 
S, so it must overlap at least one member of C and thus Z{x, y) E C, a. 
contradiction since x ^ S. So the support of an overlap class is a strong 
homogenous module. 

(2) Let A be an atom of a given overlap class C of Z{H). If A is included 
in at least two members of C, then A is exactly the intersection of all 
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members of C which include A. Since the family of homogenous modules 
is closed under intersection of overlapping members (Proposition 3), A is 
a homogenous module. Notice that if A is included in only one member 
of C, it may fail to be a homogenous module. Let us suppose that A is 
a homogenous module, and that it is overlapped by another homogenous 
module. Then it is overlapped by a homogenous module B ^ 2{H). Let 
X be an clement of B\A and y an element oi A\B. Z[x, y) contains B 
but not {y} and thus overlaps A, so it overlaps all elements of C which 
include A and thus Z{x,y) € C, a contradiction since no atom may be 
overlapped by a member of the overlap class. So the atoms of an overlap 
class which are homogenous modules are strong. 

Now, let us prove that if M is a non-trivial strong homogenous module then 
it is the support or an atom of some overlap class. We shall distinguish three 
cases. Let Mp be the strong parent of M (which exists since M ^ X). 

(1) M is prime and Mp is prime. Then for all x e M and all y G Mp \ M, 
M — Z{x,y). As M is a strong homogenous module, it alone forms a 
trivial overlap class of Z{H) and is equal to its support and to its unique 
atom. 

(2) M is prime and Mp is brittle. Then for all x e M and all y e Mp\M, M is 
included in Z{x, y). Notice that these Z{x, y) belong all to a same overlap 
class C of Z{H). Since M is a strong homogenous module of if , M C S{C) 
cannot overlap any member of C. Moreover, for all M C C S'(C), 
would overlap Z{x,y) with x E M and y E N \ M. Hence, M is by 
definition an atom of C. 

(3) M is brittle. It is easy to notice that M has k > 3 strong children 
Ml, . . . , Mfc. Let us pick an element Xi in each Mj. Then for all i and 
j we consider Z{xi,Xj). Not all of them are strong homogeneous mod- 
ules (otherwise, M would be prime). Let us consider the overlap graph of 
these homogeneous modules (the vertices are the homogeneous modules, 
and there is an edge between overlapping homogeneous modules). Each 
connected component of this graph is an overlap class. According to the 
first sentence of the theorem, the support of each overlap class is a strong 
homogeneous module. If there are two overlap classes, the support of at 
least one is a strong homogeneous module that is strictly between M and 
its sons Mj in the inclusion tree, since the overlap graph has at least one 
edge, a contradiction. So there must be only one overlap class, whose 
support is exactly M. 

□ 

For an arbitrary homogeneous relation. Theorem 2 gives the basis for an 
0(|X|"^) time enumeration of all strong homogeneous modules, which is de- 
picted in Section 4.5. 
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2-4 Particular Homogeneous Relations 



We now survey some classes of homogeneous relations defined by added ax- 
ioms, which, in practice, frequently occurs. For instance, the class of standard 
homogeneous relations (see Definition 3) has very specific properties, leading 
to efficient decomposition algorithms (see Section 5). 

Definition 5 A homogeneous relation H is said to he 

• weakly graphic if H{y\xz) A H{z\xy) =^ H{x\yz) for all x,y, z & X ; 

• weakly digraphic if H{s\xy) A H{t\xy) A H[y\sx) A H(y\tx) ^ H(x\st) 
for all X, y,s,t & X ; 

• modular quotient if H{x\st) ^ H{y\st) for all homogeneous modules 
M of H, for all x,y e M, and s,t ^ M. 

Proposition 6 A weakly graphic homogeneous relation is weakly digraphic. 
There are weakly graphic homogeneous relations that are not modular quo- 
tient. There are modular quotient homogeneous relations that are not weakly 
digraphic, hence not weakly graphic. 

Proof: If H is weakly graphic, H{s\xy) and H{y\sx) imply H{x\sy). Like- 
wise, H{t\xy) and H{y\tx) imply H{x\ty). Then, H{x\st) by transitivity of 
Hx- Hence, H is weakly digraphic. Besides, let K be defined over Xk — 
{x,y,s,t} as = {M, = {{^}As,t}}, = {{x,y},{t}}, 

and Kt = {{x,y},{s}}. Then, K is weakly graphic (exhaustive checking on 
all triplets) but not modular quotient {-iH{x\st) and H{y\st) for the homo- 
geneous module {x,y}). Finally, let L be defined over Xl = {x,y, s,t, z} 
as = {{s},{t,y,z}}, Ky = {{s,t,x}, {z}}, Ks = {{x,y}, {t, z}}, Kt = 
{{x, y}, {s, z}}. and — {{x}, {s, t, y}}. Then, L vacuously is modular quo- 
tient as having no homogeneous module, but not weakly digraphic (x, y, s, t 
form a counterexample). □ 

The modular quotient property plays an important role in modular decomposi- 
tion algorithmics. Indeed, if H is modular quotient, elements in a homogeneous 
module M of H uniformly perceive a set A not intersecting M: if one element 
of M distinguishes A then so do all. This, combined with the definition of a 
homogeneous module, allows to shrink M into a single element, the quotient 
by M, or to pick a representative element from the homogeneous module. Re- 
cursion can therefore be used when dealing with homogeneous modules. The 
modular quotient and restriction (Proposition 2) properties were first used 
in modular decomposition of graphs and are useful for algorithmics [29]. In 
this paper, these relations will be qualified as good homogeneous relations, and 
Section 5 is devoted to their study. 
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Let the congruence w.r.t. H of an element x E X stand for the number 
of equivalence classes of the relation H^. Then, the local congruence of H 
is the maximum congruence of all elements of X. Homogeneous relations of 
congruence 2 plays a special role in graph theory as they include the class of 
standard homogeneous relations of undirected graphs and tournaments (see 
next section). Furthermore, those relations satisfy the following nice property. 

Proposition 7 Any weakly graphic homogeneous relation H of local congru- 
ence 2 is modular quotient. 

Proof: Suppose H weakly graphic and not modular quotient. Then, there 
exist x, s, t pairwise distinct elements such that {x, y} is a homogeneous 
module, H{x\st), and Let us prove that we have both 

and -'H{y\xt). Indeed, suppose w.l.o.g. that H{y\xs). Then, the transitivity of 
Hy implies -iH(y\xt) (for we already have -iH(y\st)). Besides, since {x,y} is a 
homogeneous module, H{s\xy). The weakly graphic property implies H{x\sy), 
and the transitivity of yields H{x\ty). But then we would have H{x\ty), 
H{t\xy) ({x, y} homogeneous module), and -iif(|/|a;t), which is a contradiction 
with being weakly graphic. Hence, -iH{y\st), -^H{y\xs), -^H{y\xt), and the 
congruence of y is at least 3. □ 

2.5 Standard Homogeneous Relations 

Given a (directed) graph, and more generally a 2— structure, the associated 
standard homogeneous relation is defined in Definition 3. Such relations are 
peculiar and satisfy the following fundamental property. 

Proposition 8 The standard homogeneous relation of a 2— structure is mod- 
ular quotient. In particular, this result holds for graphs, tournaments, oriented 
graphs, and directed graphs. 

Proposition 8 has important algorithmic implications that will be detailed in 
Section 5. Now, the name of weakly graphic and weakly digraphic homogeneous 
relations used in the previous section is motivated by Proposition 9 below. A 
symmetric 2— structure refers to an edge-coloured clique (the clique is seen as 
an undirected graph, see e.g. [16] for further information). 

Proposition 9 The standard homogeneous relation of a directed graph, resp. 
a 2— structure, is weakly digraphic. The standard homogeneous relation of an 
undirected graph, resp. a symmetric 2— structure, is weakly graphic. 

We now investigate a converse question: given a homogeneous relation H over a 
finite set X, does there exist an undirected graph, or a tournament, admitting 
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H as standard homogeneous relation? H is defined as a graphic homogeneous 
relation if its local congruence is at most 2 and if H[{a^ b, c}] has exactly or 
2 elements of congruence 2 for every triple {a,b,c}. H is tournamental if 
its local congruence is at most 2 and if H[{a, b, c}] has exactly 1 or 3 elements 
of congruence 2 for every triple {a,b,c}. 

Theorem 3 H is the standard homogeneous relation of an undirected graph 
if and only if it is graphic. H is the standard homogeneous relation of a tour- 
nament if and only if it is tournamental. 

Proof: It is straightforward to check that the standard homogeneous relation 
of any graph, resp. tournament, is graphic, resp. tournamental. The converse 
for graphs can be proved as follows. Let if be a graphic homogeneous relation 
over a finite set X, and x & X. Let be one of the possibly two equivalence 
classes of (there always is at least one such class) . We define the matrix M 
as: M{x, y) = 1 if y E and M{x, y) = otherwise; for all x' ^ x, M{x', y) ~ 
1 if y G C^i and A'/(x, y) = otherwise, where C^' is the equivalence class of H^i 
containing x. Suppose M not symmetric. Then, there exists y ^ z both distinct 
to X such that M{y, z) ~1 and M{z, y) = 0. But then H[{x, y, z}] would have 
exactly 1 or 3 elements of congruence 2. Therefore, M is a {0, 1} symmetric 
matrix and can be seen as the adjacency matrix of some undirected graph G. It 
is then straightforward to verify that H is the standard homogeneous relation 
of G. The proof for tournaments is similar. We use the characterisation that 
the adjacency matrix of a tournament is a {—1, 1} anti-symmetric matrix since 
there are no non-edges and no double arcs. □ 

Corollary 1 It can he tested in 0{\X\^) time if a homogeneous relation H 
admits a graph G or a tournament T such that H{G) = H or H{T) = H. 

Proof: First check if all element x has congruence at most 2. Then check for 
all triples the corresponding property of the restricted relation. □ 

Notice that, if a graphic, resp. tournamental, relation H is given as \X\ sets 
of equivalence classes of (cf Section 4.1), then, the adjacency list represen- 
tation of the corresponding graph, resp. tournament, can be built in 0(|Xp) 
time. Indeed, for graphs one just has to decide which class of the first vertex 
V E X represents its neighbourhood. Then, for any other vertex u, the class 
containing v will be its neighbourhood if it is a neighbour of v, and its non- 
neighbourhood otherwise. Simply remove the "non-neighbourhood" classes (in 
0(|X|) time each): the other class in each case is the vertex's adjacency list. 
A similar construction can be performed for tournaments in the same 0(|Xp) 
worst case time. 

Remark: Extending Theorem 3 to symmetric 2— structures is quite straight- 
forward. It would be interesting to characterise the standard homogeneous 
relations of directed graphs, and 2— structures. 
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3 Pcirtitivity and Decomposition Theorem 



A generalisation of modular decomposition, known from [9] , less general than 
homogeneous relations but more powerful, is the partitive families. The sym- 
metric difference of two sets A and B, denoted by AAB, is {A\ B) Li {B \ A) . 

Definition 6 A family 3" C J'(X) is weakly partitive if it contains X and 
the singletons {x} for all x E X, and is closed under union, intersection and 

difference of overlapping members, i. e. 

A £3" A B £3" A A^B =^ AnB eS" A AUB eS" A A\B eS". 
Furthermore a weakly partitive family 5" is partitive if it is also closed under 
symmetric difference of overlapping members: 

AeS" A B eS" A A(^B =^ AAB e 3". 

Let 3" be a weakly partitive family over X. As mentioned before, strong mem- 
bers of 5" can be ordered by inclusion into a tree, so-called generalised decom- 
position tree (see Section 2.3). In this tree, the child, under the usual parental 
notion in trees, of an internal node M is by definition a strong member of 
J', which is also a strong child of the strong member M G 5", in the sense of 
Section 2.3. Besides, a weak child of the node M will refer to the definition of 
Section 2.3. Let us define three types of strong members of 5", namely three 
types of nodes of the tree: 

• prime nodes which have no weak children, 

• degenerate nodes: any union of strong children of the node belongs to 3^, 

• linear nodes: there is an ordering of the strong children of the node such 
that a union of them belongs to IF if and only if they follow consecutively 
in this ordering. 

Theorem 4 [9j In a partitive family, there are only prime and degenerate 
nodes. In a weakly partitive family, there are only prime, degenerate, and linear 
nodes. 

The generalised decomposition tree hence is an space coding of the 

family: it is sufficient to type the nodes into complete, linear or prime, and 
to order the children of the linear nodes. It is then called the decomposition 
tree of the family. From this tree, all weak members of 3" can be outputted by 
making simple combinations of the strong children of brittle (degenerate or 
linear) nodes. Now, the following property states that homogeneous modules 
of some homogeneous relations are proper generalisations of (weakly) partitive 
families. 

Proposition 10 The homogeneous modules of a weakly graphic, resp. weakly 
digraphic, homogeneous relation H form a partitive, resp. weakly partitive, 
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family. 



Proof: Proposition 3 gives the closure by intersection and union of overlap- 
ping members. Let A e M-h and B e Mh be two overlapping homogeneous 
modules of H. Suppose that there is a splitter s of A\B: there are x,y G A\B 
such that -'H(s\xy). Moreover, s E Ad B otherwise it would be a splitter of 
A. Finally, since A^^B, there exists an clement t E B \A. We have: H{x\st) 
and H{y\st) and H{t\sx) and H{t\sy) and H(t\xy). In other words, H is not 
weakly digraphic. Hence, the family of homogeneous modules of a weakly di- 
graphic homogeneous relation is weakly partitive. Besides, suppose that z is a 
splitter of AAB. Then, z E AnB and there exists x E A and y E B such that 
-iH{s\xy). Since H{x\yz) and H{y\xz), H is not weakly graphic. Hence, the 
family of homogeneous modules of a weakly graphic homogeneous relation is 
partitive. □ 

As a result, the homogeneous modules of a standard homogeneous relation 
form a weakly partitive family because such a relation always is weakly di- 
graphic (cf Section 2.5). More generally, we will prove in Proposition 15 that 
the homogeneous modules of any homogeneous relation that satisfies the mod- 
ular quotient property (cf Section 2.4), so-called good homogeneous relation, 
form a weakly partitive family. Recall that a weakly digraphic homogeneous 
relation is not necessarily modular quotient (cf Proposition 6) . 



4 Algorithms for Arbitrciry Homogeneous Relations 

This section considers a given homogeneous relation H over a ground set X, 
and builds tools for computing the generalised modular decomposition tree of 
H. The best performance to compute this tree in the general case will be given 
in 0(|X|'^) time in Section 4.5. Notice that the decomposition Theorem 4 does 
not necessarily hold in this section. 



4.1 Data Structures 

According to Definition 1, a homogeneous relation H can be represented in 
0(|Xp) space by an n x n matrix A of values in |l,n| as follows. U X = 
{xi, . . . ,Xn}, each equivalence class of the relation ifj,. will be assigned a 
distinct number from 1 to n. Then, the cell A^j has value k if and only if Xj 
belongs to the equivalence class of H^- having the value k. This representation 
allows to test in 0(1) time whether H{xi\xpXq) by checking if = Ai^g. 
However, retrieving an equivalence class requires an 0(|X|) worst case time. 
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Another alternative is to use the hst representation: each element x & X 
will be associated to a list of equivalence classes of the relation H^. This list 
is allowed to among the equivalence classes of H^, for 

instance the largest one. Thus, the total used space is 0{n + m), with n = \X\ 
and m = J2x&x{''^~ \Cx\)- Though this representation allows access in 0(1) to 
an equivalence class of for any element x, testing if H{x\yz) would require 
0{n-\C,\). 

Notice that for a homogeneous relation, it is straightforward to construct 
in 0(|Xp) time a list representation given any matrix representation, and 
conversely. 

N.B. Without further specification, all algorithms presented in this paper take 
matrix representations as input. 



4-2 Smallest Homogeneous Module Containing a Subset 

Let be a non-empty subset of X. As 1\Ih is closed under intersection, there is 
a unique smallest homogeneous module containing S, namely the intersection 
of all homogeneous modules containing S, denoted henceforth by SM{S). 



Algorithm 1: Smallest homogeneous module containing S 
Let X be an element of 5", M := {x} and F :— S \ {x} 
while F is not empty do 

pick an element y in F ; F := F \ {y} ; M :— M U {y} 

for every element z ^ (M U F) do 
\_ if ^H{z\xy) then F := F U {z} 

output M (now equals to SM{S)) 



Theorem 5 Algorithm 1 computes SM{S) m 0{\X\.\SM{S)\) = 0{\X\^) 
time. 

Proof: Time complexity is obvious as the while loop runs \M\ — 1 times 
and the for loop |X| times. The algorithm maintains the invariant that every 
splitter of M is in F. When M is replaced by MU{y}, using transitivity of the 
relation H^, every splitter for MU{y} either distinguishes x from y, or already 
is in F. The algorithm ends therefore on a homogeneous module that contains 
S, and thus we have SM{S) C M.UM ^ SM{S) let s be the first element of 
M\SM{S) added to F (eventually added to M). It distinguished two elements 
X and y from SM{S), contradicting its homogeneity. So SM{S) — M. □ 
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4-3 Maximal Homogeneous Modules Excluding an Element 



Proposition 11 Let x be an element of X . As "Mh is closed under union of 
intersecting subsets, there is a unique partition of X\{x} into Si, . . . , Sk such 
that every Si is a homogeneous module of H and is maximal w.r.t. inclusion 
in M-H- 

We call MaxM{x) this partition of maximal homogeneous modules excluding 
X, and propose a partition refining algorithm for its computation. It is straight 
from definition that 

Lemma 1 Every homogeneous module excluding x (especially the maximal 
ones) is included in some equivalence class of H^. 

Therefore our algorithm starts with the partition P = {H^, .... H^} of equiv- 
alence classes of H^. Then the partition is refined (parts are split) using the 
following rule. Let y be an element, called the pivot, and Y the part of P 
containing y. 

Rule 1 split every part A of P, except for Y, into A n Hy,. . . ,A n Hy 

Notice that a part is broken if and only if its splitters include y. 

Lemma 2 Starting from the partition Pq = {Hx--^x}> application of 
Rule 1 (for any pivot in any order) until no part can be actually split, produces 
MaxM{x). 

Proof: The refining process ends when no pivot can split a part, i.e when 
every part is a homogeneous module. Let us suppose one of these homogeneous 
modules M is not maximal w.r.t. inclusion: it is included in a homogeneous 
module M' , itself included in an equivalence class of Hx- Let us consider the 
pivot y that first broke M' . It cannot be out of M' , as M' is homogeneous 
module, nor within M' , as a pivot does not break its own part. But M' was 
broken, contradiction. □ 

Let us now implement this lemma into an efficient algorithm. Let Pj be the 
partition after the ith application of Rule 1, y be a given vertex used as pivot, 
and Yi the part of Pi containing y. We say that a part B of Pj descends from 
a part A oi Piii i < j and A G B. Clearly, after y is chosen as pivot at step 
i, y does not distinguish any part of Pj excepted 1^. If y is chosen as pivot 
after, at step j > i, y may only split the parts of Pj_i that descend from Yi. 
Only these parts have to be examined for implementing Rule 1. But Yj itself 
has not to be examined. 
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Let us suppose that, for a part A, we can split it in (9(|74|) time when applying 
Rule 1 with pivot y. Then the time spent at step j is — \Yj\), the sum 

of the size of the parts that descend from Yi save Yj. The time of all splittings 
with y as pivot is leading to an time complexity. This is 

implemented in Algorithm 2. 



Algorithm 2: Maximal Homogeneous Modules excluding x 

for every group G do 

for every part C of G do 

Compute the set Z of elements in G but not in C 
for every element y of C do 

Partition Z according to the equivalence classes of Hy 
|_ Add each partition set to the refining set pool 

Set the group boundaries to the parts boundaries (from Pj_i to Pj) 
for each refining set R of the pool do 
|_ Remove R from the pool and then refine Pj using R 



Let us suppose that the parts are implemented as a linked list [23], and the 
new parts created after splitting an old one replace it and follow consecutively 
in the list. Then for each pivot y two pointers, one on the first part that 
descends from Yi and the second to the last part, are enough to tell the parts 
to be examined. A simple sweep between the pointers, omitting Yj, gives them. 
We call all classes descending from a previous one a group. 

Now let us show how a part A can be split in 0(|A|) time. It is a classical trick 
of partition refining [22,23,31]. If the equivalence classes of Hy are numbered 
from 1 to k, then A can be bucket sorted in 0(|A| + /c) time, then each bucket 
gives a new part that descends from A. If \A\ < k, we have to renumber 
the used equivalence class of Hy from 1 to k' < \A\ before bucket sorting. A 
first sweep on A marks the used equivalence class numbers. A second sweep 
unmarks an used number the first time it is seen, and replaces it by the 
new number (an incremented counter) which is less than \A\. The vector of 
equivalence class numbers is initialised once in 0{k) time. 

The last point is the ordering in which pivots are taken. Using all elements as 
pivots, and repeating this |X| times, i.e. |Xp applications of Rule 1, is enough. 
A clever choice is to use y only if Y^ has been split, keeping a queue of "active" 
pivots. Let us define a measure that will be used later for complexity analysis. 

Definition 7 Let ^ he a partition of X . Q(7) he the number of pairs 
such that X and y are not in the same part of 7. 
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Qi?) is between 1 (for the trivial partition {X}) and (for the trivial 

partition into singletons). 

Theorem 6 MaxM{x) can be computed in Q{Q{MaxM{x))) = 0{\X\'^) 
time. 

Proof: For the correctness proof, one just has to check that the above algo- 
rithm implements correctly Lemma 2. For time complexity issues, notice that, 
for each pivot y, an element z is placed in Z only once. But it is placed in Z 
only if y and z are not in the same part. At each step, refining Z according to 
the equivalence classes of Hy, and then refining using all sets generated by y, 
takes 0(|Z|) time. Hence the algorithm takes Q{Q{MaxM{x)) time. □ 

4-4 Modular Primality test 

We recall that H is modular prime if all its homogeneous modules are trivial 
(see Section 2). 

Theorem 7 One can test in 0{\X\'^) time if H is modular prime. 

Proof: If |X| < 2 the answer is yes. Otherwise let x and y be two elements 
of X. In 0(|Xp) time, the algorithm of Section 4.3 can output the maximal 
homogeneous modules excluding x. If one of them is non-trivial then the an- 
swer is no. Otherwise all non-trivial homogeneous modules will contain x. In 
0(|Xp) time, the algorithm of Section 4.3 can output the maximal homoge- 
neous modules excluding y. If one of them is non-trivial then the answer is no. 
Otherwise all non-trivial homogeneous modules will contain x and y. Then, 
Algorithm 1 can be used with S — {x, y}, in 0(|Xp) time. The answer is yes 
if and only if SM{{x, y}) ^X. □ 

4.5 Strong Homogeneous Modules Enumeration 
Theorem 2 straightforwardly leads to an algorithm: 

Theorem 8 The strong homogeneous modules of a homogeneous relation H 
on X can he enumerated in 0{\X\^) time. 

Proof: First compute MaxM{x) for all x e X. All these sets together form 
exactly the family 2{H) defined in Theorem 2. It can be done in 0(|X|^) time 
using the algorithm of Section 4.3 |X| times. The size of this family (sum of 
the cardinals of every subset) is 0(|Xp) since they form |X| partitions. Using 
Dahlhaus's algorithm [12] the overlap components can be found in time linear 
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on the size of the family, namely According to Proposition 5 there 

are at most \X\ non-trivial overlap classes. 

For each class it is easy to compute its support, and in time easy to 

compute all its atoms. For instance, consider the vector of parts of the overlap 
class containing a given element: the atoms are the elements with the same 
vector. Sorting the list of elements of the supports 0(|X|) times, one time per 
part, gives the elements with the same vector, thus the atoms. 

Then the supports and atoms must be sorted by inclusion order 

into the inclusion tree of the strong homogeneous modules. It can be done 
in 0(|X|^) time using the same sorting technique. 

Eventually, "bad" atoms - those that are not strong homogeneous modules - 
must be removed from the tree. According to the first statement of Theorem 2, 
the atoms which are homogeneous modules are strong. We just have to perform 
0(|X|) tests on all nodes of the tree to test which of them are homogeneous 
modules, which can be done in time for each. □ 

4-6 Computation of the Generalised Decomposition Tree given a Factoring 
Permutation 

The notion of a factoring permutation in the case of graphs [7] was introduced 

to give an alternative for computing the modular decomposition tree of a 
graph without the precomputing of maximal modules excluding some vertex 
X [6,8,21,23]. It can be extended to homogeneous relation as follows. 

Definition 8 (Factoring Permutation) A factoring permutation of a ho- 
mogeneous relation refers to a depth- first search's visit order of the leaves of 
the generalised decomposition tree of the relation. 

We here address the problem of, given a homogeneous relation H over a finite 

set X and a factoring permutation a, computing the generalised modular 
decomposition tree of H . Of course the algorithm of Section 4.5 answers to this 
question. However, this section will depict a more efficient solution, 
which relates to Uno and Yagiura's iterative idea [6,33]. 

Actually, the name of factoring permutations is mainly motivated by the fol- 
lowing characterisation. Without loss of generality, we denote the elements of 
XbyX = {l,2,...,n}. 

Proposition 12 If a is a factoring permutation of a homogeneous relation H 
over a finite set X , then every strong homogeneous module of H is an interval 
of a, namely it is of the form {cr(i), a{i -|- 1), . . . , cr(j)}. 



19 



Roughly, to enumerate the strong homogeneous modules of if, it suffices to 

find among the intervals of a those that are strong homogeneous modules. Let 
lij denote the cr— interval Uj — c(i + 1), . . . , and Sy the splitter 

set of lij. 

Proposition 13 %ij = S(i+i)j U ^m+i) \ {cr(i)}. 

Proof: that S(i+i)j U S,(j+i) \ {o'(0} ^ S;? is straight from definition of a 
splitter. Conversely, let x ^ lij be such that x ^ S(i+i)j U Si(i+i). Then, by the 
transitivity property of H^, we obtain H{x\yz) for all y,z & lij, or in other 
words X ^ Sjj. Hence, §ij C S(i+i)j U Si(i+i). We use the fact that a{i) ^ Sjj to 
conclude. □ 

This leads to a naive 0(|X|^) solution to this section's question: for all in- 
terval lij, compute Si(i+i), then $ij using the previously computed §(i+i)j and 
Proposition 13, eventually test if Sij is empty. Let us now improve this idea. 
The interval lij is said to be right-free if it does not have a splitter on the right 
in the order a, namely for all k > j, a{k) does not belong to $ij. Obviously, if 
lij is a strong homogeneous module, lij is right-free. However, a much more 
interesting viewpoint is as follows. If lij is not right-free, then there will be no 
i' <i such that liij is a strong homogeneous module. Furthermore, 

Proposition 14, If ji < ■ ■ ■ < are such that any lij^ {1 < q < k) is 
right-free, then $ij^ C • • • C §ij^ . 

Proof: All splitters of these intervals stand on the left of a{i) in the order a. 
Hence, a splitter s of lij^ can not belong to /ij^+i , and will belong to Sij^+i . □ 

Roughly, if in some iteration step 1 < i < n, we only store some right-free 
intervals in a list RF = {hj^, . . . jlij^), then all their corresponding splitters 
can easily be stored by differences in a list AS* = (A^j, . . . , Aj^), where Aj^ = 
Sjjj and Aj^ = Sj^^ \ §ij^_^ > 2). Under this convention, an interval lij^ of 
the collection is a homogeneous module if and only if all the q*^ first members 
of A^" are empty: A^-^ = • • • = A^-^ = 0. 

From iteration step i to {i — 1), the collection of intervals will extend from 
RF = {lij^, lij^) to RF = . . . , /(i-i)jj, and the list AS 

will be updated accordingly using Proposition 13. Also, if for some jg, the 
extension of lij^ to I{i-i)jq introduces a splitter a{k) such that k > jg, then we 
remove this interval from RF for it no more is right-free and jg will have no 
chance to be the right boundary of an unvisited strong homogeneous module. 
We come to Algorithm 3. For convenience, each interval lij^ will be represented 
by its right boundary: we shall use RF — {ji, ■ ■ ■ ,3k)- 

Remark: Basically, the first step i = noi the main loop still is an initialisation 
step: at the end of the loop, we always have RF = (n), A^ = (0), and 
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Algorithm 3: Generalised modular decomposition tree computation from a 
factoring permutation 

Input: a homogeneous relation H over a finite set X, and a factoring 

permutation a of H 
Output: the generalised modular decomposition tree T oi H 
RF ^ and ^ () and M ^ 

Create a dummy y = a{n + 1) such that H{s\xy) for all s e X and x — a{n) 
for i — n downto 1 do 

X ^ a{i) and y ^ a{i + 1) 

if X belongs to some member of AS then remove x from that member 
for every s = a{l) with I < i and -'H{s\xy) do 
L Add s to the first member of AS 

Find s = ij{r) such that -iH{s\xy) and r maximum, otherwise r 
while the first member j of RF satisfies j < r do 
Remove j from RF 

Let Fst and Snd be the first and second members of AS 
L Snd ^ Snd U Fst and remove Fst from AS 
RF ^ {i, RF) and AS ^ (0, AS) 
Let S, resp. j, be the first member of A^, resp. RF 
while 5" = do 
M ^ {lij} U M 
|_ Let S, resp. j, be its next member in AS, resp. RF 

Remove the weak members of M 

Construct T, the inclusion order of members of M 

Output T 



M — {n}. The real computation starts at step i — n — 1. 

Inv£triant 1 For alll < i < n, let RFi = (ji, . . . , jk) and ASi = (Ai, . . . , A^) 
be the values of RF and AS, at the end of the first loop "for" in Algorithm 3. 
Then, 

• for all member j of RFi, the interval lij is right-free; 

• for alll<q<k, = Ai U ■ ■ ■ U Ag. 

Algorithm 3 correctness directly follows from Invariant 1. As for complexity 
issues, it is quite straightforward to check that the computing time of all loops 
is in 0{n^). After those loops, removing weak members of the list M can be 
done in linear time on \ M\ using the lexical member ordering of M: lij is before 
li/ji in M if and only if i < i' or (z = z') A (j < /). Notice that \M\ is less than 
the number of intervals of a, which is in 0{n^). Likewise, the time spent for 
ordering by inclusion the remaining members of M is linear on their number 
using the lexical property. Whence, the global computing time of Algorithm 3 
is 0{n^). 
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Theorem 9 Given a factoring permutation a of a homogeneous relation H 
over a finite set X , one can compute the generalised modular decomposition 
tree of H in 0(|Xp) time. 

Factoring permutations can be get in 0(|Xp) time in many cases, especially 
with standard homogeneous relations of 

• inheritance graphs: a linear extension gives a factoring permutation [14]; 

• chordal graphs: the cardinality lexicographic breadth first search of the 
graph yields a factoring permutation [24] ; 

• tournaments: a very simple partition refining algorithm (greedily choose x 
and partition the class containing x into N~{x), {x}, N~^{x)) computes a 
factoring permutation [26]; 

• undirected graphs: more sophisticated algorithms run in O(mlogn) time 
[23] or 0{n + m) time [21]. 



5 Good Homogeneous Relation Decomposition Algorithm 

The good homogeneous relations refer to homogeneous relations fulfilling the 
modular quotient property (cf Section 2.4). For instance, standard homoge- 
neous relations are good (Proposition 8). Their study is motivated by, among 
others, the following essential property. 

Proposition 15 The homogeneous modules of a good homogeneous relation 
form a weakly partitive family. 

Proof: Proposition 3 gives the closure under intersection and union of over- 
lapping members. We just have to check that, for two homogeneous modules 
A and B of H, if A<^B then A\B is a. homogeneous module. Let us suppose 
A \ B has a sphtter s. As A is a homogeneous module, s E A (1 B. Let x 
and y be two elements of A \ B such that -iH{s\xy). As A'^B there exists 
t E B \A. Since B is a homogeneous module, the modular quotient property 
gives -ii?(t|a;y). But then A no more is a homogeneous module. □ 

Let if be a good homogeneous relation over a finite set X. We address the 
problem of computing the modular decomposition tree of H, namely the in- 
clusion order of strong homogeneous modules of H. Here again, the algorithm 
of Section 4.5 can be used to give a solution to this question in 0(|A|^) time. 
However, this section will give a more efficient 0(|A|^) time solution, which 
is inspired from Ehrenfeucht et al. works [15]. 

Definition 9 A super-modular-decomposition-tree (SMDT for short) of a 
good relation H on X is a tree 
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• where the leaf-set is X 

• such that each node of the tree is a homogeneous module of H 

• such that each strong homogeneous module of H is a node of the tree. 

The idea of the algorithm is to compute the left branch ( "caterpillar" ) of a 
super modular decomposition tree of H, going from the root X to an arbitrary 
element x (see Fig. 4). Then, the algorithm recurses to compute the "legs" 
of the caterpillar, and appends them to the caterpillar. Algorithm 4 captures 
this idea. Eventually, the SMDT is cast into the modular decomposition tree. 

Proposition 16 Algorithm 4 computes a super modular decomposition tree 

Proof: Obviously all outputted nodes are homogeneous modules. We just 
have to check that the tree contains all strong homogeneous modules. This is 
true indeed, because, for a strong homogeneous module M, the first element 
a; e M taken for the branch (see the definition below) at some recursive 
step outputs M G !B(a;). The goodness of the the relation gives that, when the 
algorithm is applied recursively on H[N] and when is a homogeneous mod- 
ule, the homogeneous module M of H[N] output is exactly the homogeneous 
module M oi H. □ 

Algorithm 4: Super Modular Decomposition Tree of a Good Homogeneous 

Relation 

Input: a good homogeneous relation H over a finite set X 
Output: a super modular decomposition tree T oi H 
Let X be an element of X 

Compute MaxM{x), the maximal homogeneous modules not containing x 

Order MaxM{x) = Mi..Mk such that for each Bj e'B, 1 <j <l, there 

exists f{j) such that Bj = {a;} tb) Mi tt) M2... tt) M/(j) 

Initiahse T to be the a;— branch 

for every (1 < i < k) do 

Compute recursively the modular decomposition tree % of H[Mi\ 
Append % to the node Bj of T such that j < i < f{j) 

Output T 

We are now to give a solution to each step of Algorithm 4, and prove their 
correctness. 



Fig. 4. A recursive approach to compute a super modular decomposition tree. 



x-branch 




recurse here 
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5. 1 Strong homogeneous modules containing x 



Definition 10 (x— branch) The x— branch of a good homogeneous relation 
H over X is the set of all strong homogeneous modules containing the 

element x & X , ordered by inclusion. In other words, it is the path from the 
root to leaf x of the modular decomposition tree of the relation. 

The tool to construct the strong homogeneous modules containing x is the 
construction of the maximal homogeneous modules excluding x. Section 4.3 
defined the set MaxM{x) = {Mi, . . . , M^} of maximal homogcncoTis modules 
excluding x, which is a partition of X \ {x} by Proposition 11. Let us examine 
the relationship between MaxM{x) and 'B{x) 

Proposition 17 The homogeneous modules of MaxM{x) — {Mi,...,Mk} 
can be ordered from 1 to k in such a way that 

for each B e 'B{x), there exists f such that S = {x} W Mi tt) M2... tb) Mf. 

Proof: For a homogeneous module B G 23 (x), the maximal homogeneous 
modules not containing B form a partition of X. Of course each homogeneous 
module of MaxM{x) is included (or equal to) one of the homogeneous modules 
of this partition. So a homogeneous module of MaxM{x) can not overlap 
a homogeneous module B e ^(a;), and the proposition follows. Indeed, to 
construct the ordering, just number the homogeneous modules of 'B{x) from 
Bq = {x} to Bi = X using inclusion order. Then number the homogeneous 
modules of MaxM{x) included in Bi from 1 to /(I), the homogeneous modules 
of MaxM{x) included in B2 but not in Bi from /(1) + 1 to /(2), and generally 
the homogeneous modules included in B^ but not in from f{i — 1) + 1 to 

m- □ 

A consequence is that, if we order the elements of the a;— branch from Bq = {x} 
to Bi — X in increasing inclusion order, then for all 1 < i < Z Sj+i \ Bj, is 
equal to some elements of MaxM{x) that follow consecutively in the above 
ordering. The following fact is obvious. 

Proposition 18 Let B^ e 'S>{x) be a non-leaf strong homogeneous module 
containing x, Cl...Cf^^^ be its children in the modular decomposition tree and 
j such that Cf = -Bj+i is the child containing x. If Bi is linear we suppose the 
children are ordered according to the linear ordering. 

• If Bi is prime then for all k ^ j Cf e MaxM{x) 

• IfBi is linear then {jlZt^ e MaxM{x) and Utf+i C'f e MaxM(x) 

• If Bi is complete then [JkytjCi G MaxM{x) 

There are no more elements in MaxM{x) than those described above. 
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5.2 Quotient relation 



Now let us construct a quotient relation. For all G MaxM{x) let Cj e Mj 
be a representative element of Mj (an arbitrary element). The quotient relation 
of H by MaxM{x), denoted H{x), is the relation 

H{x) = H[{x,ei, . . .,ek}]. 



Proposition 19 The quotient relation of H by MaxM{x) does not depend 
on the choice of the representative elements for each M^. 

Proof: This is because the relation H is good. □ 

Proposition 20 Every non-trivial homogeneous module of H{x) contains x. 

Proof: Suppose there is a non-trivial homogeneous module Ujg/{ej} of H{x) 
that excludes x. Then, |/| > 2, and Ui^jMi is a homogeneous module of H 
that excludes x, larger than an element of MaxM{x), a contradiction. □ 

For Mi G MaxM{x), let S{Mi) G 23 (x) be the smallest homogeneous module 
of !B(a;) containing Mj. Using the notations of Proposition 17 if S{Mi) = Bj 
then i < j < Proposition 18 gives the relationship between Mj and S{Mi) 
with respect to S{Mi) type (complete, linear or prime). We say that G Mj 
is a P-element (resp. L-element, C-element) if S{Mi) is prime (resp. linear, 
complete). Two elements Ci G Mi and cj G Mj are companion one of each 
other if S{Mi) = S{Mj). Proposition 18 tells that Cj has zero companion if 
S{Mi) is complete, zero or one if 5'(Mj) is linear and at least one if S{Mi) is 
prime. 

5.3 Forcing graph 

Definition 11 (Forcing Graph) Keeping the above notations, the directed 
forcing graph G{x) = {V,A) is defined as V = (ci, . . . , e^}; and an arc 
{ei,ej) G A exists if and only if —>H{ej\x,ei). 

Proposition 21 Let y be a vertex of G{x) and N*{y) the descendants of y in 
G{x) (including y itself). N*{y) U {x} is the smallest homogeneous module of 
H{x) containing y. 

Proof: First notice that all nontrivial homogeneous modules of H{x) contain 
X. Then, if the forcing graph has an edge (cj, ej) then any nontrivial homoge- 
neous module of H{x) containing also contains ej. All descendants of y in 
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G{x) are thus in any homogeneous module containing y (and x). 

Now we shall prove that for any set A of vertices of G{x) with no outgoing 
arc, A U {x} is a homogeneous module of H{x). Indeed, for all u G A and all 
v^Awe have H{v\x, m). As if is a transitive relation, then for all u,u' E A 
H{v\u,u') and thus A U {x} is a homogeneous module. So N*{y) U {x} is a 
homogeneous module of H{x). □ 

Let C be a strongly connected component (SCC for short) of G{x). The above 
proposition gives that all vertices of C are companions. Furthermore we have: 



Proposition 22 A non-trivial strongly connected components of G{x) is 
formed by companion P- elements. Conversely a maximal set of companion 
P- elements is strongly connected. 

Proof: According to Proposition 18 there are no companion C-elements and 
at most two companion L-elements. But clearly there is no arc between them. 
So a SCC with at least two vertices contains companion P-elements. Accord- 
ing to Proposition 21 if companion P-elements were split into two (or more) 
SCC C and D, then there would be either a homogeneous module of H{x) 
containing C but not D, or a homogeneous module of H{x) containing D but 
not G. In both case, the smallest homogeneous module of H{x) containing 
G U D can not be prime. □ 

Proposition 23 Two companion L-elements are false twins (they share the 
same neighbourhood and there is no arc between them). Conversely the pairs 
of false twins are exactly the companion L-elements. 

Proof: Let e and e' be two companion L-elements. The smallest homogeneous 
module M of H{x) containing {e, e'} is thus a linear homogeneous module 
{e} U M' U {e'} where M' is the strong homogeneous module son of M in the 
modular decomposition tree of H(x). Of course x E M'. Both {e} U M' and 
M' U {e'} are homogeneous modules, and the descendants of e are exactly the 
descendants of e' and are M', according to Proposition 21. Furthermore since 
H is good, e and e' are twins. □ 

According to the Propositions 21, 22 and 23 we have: 

Proposition 24 Any linear extension (topological sort) of G{x) will order 
MaxM{x) into the ordering of Proposition 17. 

Proposition 25 The x— branch of H can be computed in 0{Q{MaxM{x))) 
time 

Proof: Remind that (5( ?) is the number of pairs {a;, y} whose vertices are not 
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in the same part of a partition CP (Definition 7). Let k be the number of parts of 
MaxM{x). obviously k"^ = 0{Q{MaxM{x))) and Q{MaxM{x)) = 0(1^^^). 
The algorithm is 

• The maximal homogeneous modules MaxM{x) excluding x can be com- 
puted in time 0{Q{MaxM{x))), according to Theorem 6, using the algo- 
rithm of Section 4.3. 

• Then, the vertices of the forcing graph are determined arbitrarily: for all 
1 < i < A; let e M^. 

• Then, constructing the forcing graph G{x) in 0{k'^) time is obvious 

• Then, the topological sort G{x) in 0{k'^) time is also easy. 

• Lastly Proposition 17 tells how the ordering of MaxM(x) allow to con- 
struct 'B{x). Notice that all companion vertices appear consecutively in the 
topological sort and are all regrouped to form \ Bi. 

□ 

We thus have: 

Theorem 10 Algorithm 4 computes a super homogeneous modular decompo- 
sition tree in 0{\X\'^) time. 

Proof: This is a direct application of Propositions 16 and 25. We just have 
to show that the sum of all 0{Q{MaxM{x))) time computations is 0(|Xp). 
This is true because Q{MaxM{x)) is the number of pairs {x,y} belonging 
to two elements of MaxM[x). As the algorithm is recursively launched on a 
homogeneous module of MaxM{x), each pair {x,y} is counted once, in the 
recursive call of its least common ancestor of the SMDT finally output. □ 



5.4 Testing for weak homogeneous modules and typing the nodes 

Now, by constructing recursively x— branches, we can build a super homo- 
geneous modular decomposition tree. This tree however is not the modular 
decomposition tree of H since: 

• Its nodes are not typed complete, linear or prime, 

• It contains all strong homogeneous modules but may also contain weak 
homogeneous modules. 

Definition 12 Let N be a node of a SMDT of H, with sons Si, ... , Sk, and 
Cj e Si be an arbitrary element. The quotient oi H hy N is H[{ei, e^}]. 

Proposition 26 The quotient relation of a node N of a SMDT is either 
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• type P: with no non-trivial homogeneous module, 

• type L: the elements can he linearly ordered in such a way that the homo- 
geneous modules of the quotient relations are exactly the intervals of the 
relation, 

• type C: every subset is a homogeneous module. 

If N has k sons, a trivial 0{k'^) time algorithm can test the type and order 
the elements if needed. A classical (and easy to prove) result is that 

Proposition 27 Let T be a tree with n leaves and no node with only one 
child. Then 

J2 degree{Nf = 0{n^). 

N node of T 

We can therefore perform quadratic-time computations on each node of a 
SMDT. A first apphcation of Proposition 27 is 

Proposition 28 Let H he a good relation on X. It take time to 

compute the quotient relations for all nodes of a SMDT of H. 

Proof: A bottom-up sweep, keeping one representative per child, builds the 
representatives. Each quotient relation can then be computed in time linear 
on its size, i.e. 0{k'^). □ 

A second application of Proposition 27 together with Proposition 26 gives 
that the typing of the nodes of a SMDT takes 0(|Xp) time. Note that we 
abusively consider that weak homogeneous modules have a type. Then we can 
look for the weak homogeneous modules, and cast the SMDT into the genuine 
modular decomposition tree, using: 

Proposition 29 Let H he a good relation on X , and N he a node of a SMDT, 
and F he its father in the SMDT. F has another son A. If F is linear then 
take A that immediately precedes or follow N in the linear ordering. Take an 
element a & A. If N is non-trivial it has at least two sons B and C. If N is 
linear then take B its first child and C its last child. Finally take b & B and 
c & C. Then 

N is a weak homogeneous module if and only if 

{a,b} or {a,c} is a homogeneous module of H[{a,b,c}]. 

Proof: If {a, b} or {a, c} is a homogeneous module N is obviously weak. Con- 
versely if N is weak then it is overlapped by a homogeneous module N'. N 
and N' have thus the same father F in the modular decomposition tree. If F 
is complete, any union of a son of F included in plus one not included A^' 
overlaps N. As b & N and a G (-F \ A^) we get the result. And if F is linear 
(any other arc is excluded), then either the first son of F included in A^ plus 
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the preceding one in the hnear order, overlaps N, and {a, b} is a homogeneous 
module; or the last son of F included in N plus the following one in the linear 
order, overlap N, and {a, c} is a homogeneous module. □ 

This proposition, together with a third application of Proposition 27, gives 
that the weak homogeneous modules can be removed from a SMDT in 
time. We finally have 

Proposition 30 A Super Modular Decomposition Tree of H can be cast into 
the modular decomposition tree of H in time. 

And, together with Theorem 10 we have: 

Theorem 11 The modular decomposition tree of a good relation H over X 
can be built in 0{\X\'^) time. 

Conjecture: When the homogeneous relation is given by list representation 
(see Section 4.1), the decomposition tree can be built in 0(n + mlogn) time, 
where n = |X| and m the total length of the lists in this representation. 



6 Outcomes 

Let us examine in the sequel some of the applications of this homogeneity 
theory to modular decomposition of graphs and 2-structures, and to other 
graph relations. 

From Proposition 9 and Section 3, the modules of an undirected graph and 
of a symmetric 2-structure form a partitive family, while the modules of a 
directed graph just form a weakly partitive family. All know properties of 
modular decomposition [29] can be derived from this result. An O(n^) modular 
decomposition algorithm can also be derived from Section 5 algorithm. It runs 
in optimal time for relations given as matrices (like an adjacency matrix), but 
it is less efficient than the existing algorithms for graphs stored using adjacency 
lists [6,8,11,13,15,22,23,27]. 

In a graph we can consider different homogeneous relations, for instance the 
relation "there exists a path from vertex x to vertex y avoiding the vertex s 
or a more general relation "there exists a path from x to y avoiding the neigh- 
bourhood ofs". It is easy to see that these two relations fulfil the basic axioms 
(symmetry, reflexivity and transitivity). In the first case, the strong hommo- 
geneous modules form a partition (into the 2-vertex-connected components, 
minus the articulation points). The second relation is related to decomposition 
into star cutsets. 
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Another interesting relation is Dk{s\xy) if d{s,x) < k and d{s,y) < k, where 
d{x, y) denotes the distance between x and y. The case k = 1 corresponds to 
modular decomposition. It is worth investigating the general case. 



7 Conclusion 

We hope that this homogeneity theory will have many other applications and 
will be useful to decompose automata [1] and boolean functions [3]. Obvi- 
ously, the algorithmic framework presented here can be optimised in each 
particular application, as it has been done for modular graph decomposition 
[6,8,11,13,15,22,23,27]. 
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